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Abstract 

An approximation Ansatz for the operator solution, U(z',z), of a hyper- 
bolic first-order pseudodifferential equation, d z + a(z, x, D x ) with Re(a) > 
0, is constructed as the composition of global Fourier integral operators 
with complex phases. An estimate of the operator norm in L(H^ 3 \ H^') 
of these operators is provided which allows to prove a convergence result 
for the Ansatz to U(z' , z) in some Sobolev space as the number of opera- 
tors in the composition goes to oo. 

AMS 2000 subject classification: 35L05, 35L80, 35S10, 35S30, 86A15. 

Introduction 

We consider the Cauchy problem 

(0.1) d z u + a(z, x, D x )u = 0, 0<z<Z 

(0.2) u | z=0 = it , 

with Z > and a(z,x,^) continuous with respect to (w.r.t.) z with values in 
5 1 (M™ x R"). We denote D x — ±d x . Further assumptions will be made on the 
symbol a(z, x, £). When a(z, x, £) is in fact independent of x and z it is natural 
to treat such a problem by means of Fourier transformation: 

u(z, x') = J J exp[i(x' — x\C) — za(()] uq(x) ct£ dx, 

where <f£ := (^) n d£. For this to be well denned for all u G J^(R") or some 
Sobolev space we shall impose the real part of the principal symbol of a to be 
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non-negative. When the symbol a depends on both x and z we can naively 
expect 

u(z, x') ~ Ui(z, x 1 ) exp[i(x' — x\£) — za(0, x', £)] uq{x) et£ dx 



for z small and hence approximately solve the Cauchy problem (|0.1|) - (|0.2[l for 
z £ [0, z^] with z' 1 ) small. If we want to progress in the z direction we have 
to solve the Cauchy problem 

d z u + a(z, x, D x )u = 0, z^ < z < Z 

u ( z , •) \z=zW = u l( z(1) i ■)> 

which we approximatively solve by 

u(z, x) ~ U2(z, x) 

exp[i{xf-x\$ -{z- z {1) )a{z [1 \x\S,)] Ul (z {1 \x) dx. 

Such a procedure then goes on. 

If we call G( z ',z) the operator with kernel 

G(z',z){x' ,x) = J exp[i(x' -x\£)]exp[-(z' - z)a(z,x', 

we then see that the procedure described above involves composing such oper- 
ators: with chosen values < z^ 1 -* < • • • < z^") < z we then have 

u k+1 (z,x) = G(z,zM) ° ° ••• £(z<i>,o)( u o)(aO, 

for z > z( k \ We then naturally define for = {z^°\ z^\ . . . , z^}, a subdivi- 
sion of [0, Z] with = < z^ < ■ ■ ■ < = Z, the operator W<p,z 



Q(z,o) if < z<z«, 
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G(z,zW)l[G(zW,z(*-V) if z(k) ^ Z ^ z(k+1) - 



The procedure described above yields Wsp, z (uo) as an approximation Ansatz for 
the solution to the Cauchy problem l|U.l[l - Hl).2[l . We denote Am = sup i=1 n{zi— 
The operator Q/ Z i x \ is often referred to as the thin-slab propagator (see 
e.g. 00). 

Note that a similar procedure can be used to show the existence of an evolu- 
tion system by approximating it by composition of semigroup solutions of the 
Cauchy problem with z 'frozen' in a(z,x,D x ) [191 111) . Note that the thin-slab 
propagator Gr z > lZ ) is however not a semigroup nor an evolution family here (see 
Section for simple arguments). 

The approximation Ansatz proposed here is a tool to compute approximations 
of the exact solution to the Cauchy problem (|U.1|) - (|U.2|1 . Such computations in 
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the case of geophysical problems can be found in [3] . In exploration seismology 
one is confronted with solving equations of the type 

(0.3) (8 Z - ib{z, x, D t , D x ) + c(z, x, D u D x ))v = 0, 

(0.4) w(0, .)=«+((),.), 

where t is time, z is the vertical coordinate and x is the lateral or transverse 
coordinate. The operators b and c are of first order, wiht real principal parts, b\ 
and ci, and c\ (z, x, r, £) is non-negative. We see that the Cauchy problem ({0.1(1 - 
(|0.2(l studied here is more general. The Cauchy problem 1(0. 3[1 - 1(0.4(1 is obtain 
by a (microlocal) decoupling of the up-going and down-going wavefields in the 
acoustic wave equation (see Appendix lAl and [23 f° r details). The proposed 
Ansatz can then be a tool to approximate in practice the exact solution of 
the Cauchy problem 1(0.3(1 - 1(0.4(1 for the purpose of imaging the Earth's interior 
E] ■ As explained in Appendix the operator c works as a damping term 
that suppresses singularities in the microlocal region where its symbol does not 
vanishes. This effect is actually recovered in the proposed Ansatz. Geophysists 
are not only interested in the convergence of this Ansatz to the exact solution 
of the Cauchy problem 1(0.3(1 - 1(0.4(1 but they expect the wavefront set of the 
approximate solution to be close, in some sense, to that of the exact solution 
because seismic imaging aims at imaging the singularities in the subsurface 
(see for instance [231 E|) • We shall investigate the microlocal properties of the 
proposed Ansatz in Part II, written in collaboration with Giinther Hormann. 

In the present paper, we are interested in the analysis of the convergence of the 
approximation scheme Wqj in Sobolev spaces. Section ^ introduces the Cauchy 
problem we study and the precise assumptions made on the symbol a(z, x, £), 
especially on the real part, ci, and imaginary part, — b\, of its principal symbol. 
In Sectional we shall at first concentrate our study on the operator Q/ Z i^ z \, yet 
to be properly defined. Under some assumptions on a(z, x, £), we shall prove 
that G(z',z) is a global Fourier integral operator (FIO) with complex phase and 
that it maps 5? into J^, S?' into ,9" and ifW m to for any s. An estimation 
°f ||^(z',#) ll(fl"C<0 w ill be the first step towards the analysis in Section [3] of 
the convergence of W<p iZ . In fact we prove that for z' — z sufficiently small then 
(Theorem E23 

l|£v,,)ll(HC*),i?(.)) <l + \z'-z\M, 

for some constant M. Such an estimate is achieved by the analysis of the 
behavior of the symbol exp[— Aci] as an element of S"? , in particular as A = z'— z 

2 

goes to zero. 

In Section we study the convergence of the Ansatz W<p, z (uo) to the solution 
of the Cauchy problem (|0.1|I - H0.2J) in Sobolev spaces as A<p goes to 0. A conver- 
gence in norm of W<p iZ to the solution operator of the Cauchy problem H().l(l - I|0.2(l 
is actually obtained fTheorem l3.1Ufl : 

lim \\W V , Z - U(z,0)\\ (H (s+i) >HM) = 0, 

Zi<p — >u 

with a convergence rate of order ^ . 

At the end Section |3| we relax some regularity property of the symbol a(z, .) 
w.r.t. z by the introduction of another, yet natural, Ansatz: following |17) . the 
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thin-slab propagator, G( z ',z)i is replaced by the operator G( z ',z) with kernel 
G( z > tX )(x',x) = J exp[i{x' - x\g)] exp[-f* a(s,x' ,£)ds\ <f£. 

In Part II, we shall focus on the microfocal aspects of the operator W<p, z and 
how it propagates the singularities of the initial condition uq- We shall show 
that the wavefront set of W<p ;Z (wo)(.z, .) converges in some sense to that of the 
solution u(z, .) of the Cauchy problem (|0.1fl - HU.2f> as ^tp goes to 0. 

Multi-composition of FIOs to approximate solutions of Cauchy problems where 
first proposed in JHj and ^H] where the exact solution operator of a first order 
hyperbolic system is approximated with a different Ansatz, up to a regularizing 
operator. The technique is based on the computation and the estimation of the 
phase functions and amplitudes of the FIO resulting from these multi-products, 
a result know as the Kumano-go-Taniguchi theorem. The technique was then 
further applied to Schrodinger equations with specific symbols |12U17j . In these 
latter works the multi-product in also interpreted as an iterated integral of 
Feynman's type and convergence is studied in a weak sense. In |12) a convergence 
result in L 2 is proved. This is the type of results sought here for first order 
hyperbolic equations. We however do not use the apparatus of multi-phases 
and rather focus on estimating the Sobolev regularity of each term in the multi- 
product of FIOs in the proposed Ansatz. While the resulting product is an FIO, 
we do not compute its phase and amplitude. 

In this paper, when the constant C is used, its value may change from one line 
to the other. If we want to keep track of the value of a constant we shall use 
another letter. When we shall write that a function is bounded w.r.t. z and/or 
A we shall actually mean that z is to be taken in the interval [0, Z] and A in 
some interval [0, A maa; ] unless otherwise stipulated. We shall generally write X, 
X', X", JW, . . . , XW for R n , according to variables, e.g., x, x' , ... , x {N \ 

In the present paper, symbol spaces are spaces of global symbols; a function 
a £ tf°°{W l x W) is in S™ s (R™ x MP), < p < 1, < 8 < 1, if for all 
multi-indices a, (3 there exists C a p > such that 

\d™dla{x,0\ <C Y a /3 (l + |e|) m " p|/3|+5|Q| , xel", Cer 
The best possible constants C' a p, i.e., 

Pa p(a):= sup (l + \S\)- m+ PW-<M\d2dla(x,Q\, 

define seminorms that turn S™ s (R n x W) into a Frechet space. As usual we 
write S™ (E n x W) in the case p = 1 - 8, \ < p < 1, and 5 m (R" x MP) in the 
case p = 1, 8 = 0. 

We shall use, in a standard way, the notation # for the composition of symbols 
of pseudodiffcrential operators (-0DO). When given an amplitude p(x,y,£) £ 
SJj (X x X x R") , p > 8, we shall also use the notation a {p} (x, £) for the symbol 
of the pseudodifferential operator with amplitude p. For p £ S™ S (X x R n ) 
we shall write p* for the symbol of the adjoint operator. When composing 
■0DOs or computing adjoints of ^DOs we shall make use of the oscillatory 
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integral representation of the resulting symbol instead of asymptotic series for 
two reasons. First, we aim at estimating operator norm in L(H S , H s ) while using 
asymptotic series representations for symbols yields results up to regularizing 
operators which operator norms cannot be controlled. Second, we shall consider 
symbols in S 1 ™, for some m, including the case p = \ for which the asymptotic 
formulae of the calculus of -0DOs cease to hold. 

For r £ R we let E^ be the ^DO with symbol (£) r := (1 + |£| 2 ) r/2 . The 
operator E^ maps H^(X) onto H^-'^iX) unitarily for all s £ R with E^ 
being the inverse map. 

1 The homogeneous first-order hyperbolic equa- 
tion 

Let s £ R and Z > 0. We consider the Cauchy problem 

(1.5) d z u + a(z,x, D x )u = 0, < z < Z, 

(1.6) H-o = «o£i? (s+1) (n 
where the symbol a(z, X, £) satisfies the following assumption 
Assumption 1.1. 

a z (x, = a(z, x, £) = -z 6(z, x, £) + c(z, x, £) 

where b £ ^ ([0, Z], S 1 (M n x 1")), with real principal symbol b\ homogeneous of 
degree 1 for |£| large enough and c £ ([0, Z], S 1 (W n x R ra )) with non-negative 
principal symbol c\ homogeneous of degree 1 for |£| large enough. Without loss 
of generality we can assume that b\ and C\ are homogeneous of degree 1 for 

l£l>i. 

In Section we shall further make the following assumption. 

Assumption 1.2. The symbol a(z,.) is assumed to be in ([0, Z], S 1 (M n x 
R™)), i.e. Lipschitz continuous w.r.t. z with values in S' 1 (R ra x R"), in the sense 
that, 

a{z , x, £) — a(z, x, £) = (z — z)a(z , z, x, £), < z < z < Z 
with a(z',z,x,£) bounded w.r.t. z' and z with values in S 1 (W l x R ra ). 

Weaker assumptions will also be formulated in Section OH for instance by the 
introduction of another approximating Ansatz. 

We denote by a± = —ibi + c± the principal symbol of a and write b = b\ +&o with 
b £ Sf°([0, Z],S°(R n x R")) and c = c x + c with c £ Sf°([0, Z], S°(R n x R")). 
Assumption 11.11 ensures that the hypotheses (i)-(iii) of Theorem 23.1.2 of [B] 
are satisfied. Then there exists a unique solution in ^°([0, Z], i2"( s+1 ) (R™)) n 
^([0, Z], #< S >(R")) to the Cauchy problem :fO; fPti. 
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Furthermore, we have the following energy estimate Lemma 23.1.1] for any 
function in ^([0, Z], tf( s >(R™)) n ^°([0, Z], H^+V (R n )) 

(1.7) sup exp[-Az] \\u(z,.)\\ H(B) < ||u(0,.)|| ff ( a ) 

z£[0,Z] 

+ 2 / exp[-Az] \\d z u + a z (x, D x )u\\ H{s) dz, 
Jo 

with A large enough (A solely depending on s). 

By Proposition 9.3 in [SJ Chapter VI] the family of operators (a z ) z& \o z\ gen- 
erates a strongly continuous evolution system. We denote U(z', z) the corre- 
sponding evolution system: 

U(z", z') o U(z', z) = U(z", z), Z > z" > z' > z > 0. 
Then the Cauchy problem (|1.5[1 - H1.6|) reads 

d z U(z, zq)uq + a(z, x, D x )U(z, zq)uq — 0, < zq < z < Z, 

[/(z„,z )u„ = Wo a (s+1) (i") 

while U(z, z )u e H^ s+ ^(R n ) for all z e [z , Z). 

2 The thin-slab propagator. Regularity proper- 
ties. 

We follow the terminology introduced in Sections 25.4-5] for FIOs with 
complex phase. Let z',z £ [Q,Z] with z' > z and let A := z' — z. Define 

<P(z',z) € tf°°(X' x X x M") as 

(2.8) <t) {zl , z) {x',x,() := {x' -x\0+iA ai {z,x',0 

= {x' - x\Q + A6i(«y ,£) +iA Cl (z,x' ',£). 

Remark 2.1. The function 4>ut tZ ) is assumed to be homogeneous of degree 1 
only when |£| > 1. This however is not an obstacle to the subsequent analysis, 
e.g., FIO properties, since to define such operators the phase function need not 
be homogeneous of degree 1 for small |£|. In the subsequent results concerning 
the phase function and FIOs one will then assume that |£| is large enough, i.e., 
I£l > I- 

Lemma 2.2. cj>( z > tZ ) is a non-degenerate complex phase function of positive type 
(at any point (x' ,xo,^o) where d^4>( z ', z ) = 0). 

Proof. Note that, by Assumption 11.11 Im((f>( z t z \) > and <f> is homogeneous 
of degree 1; d x <fi = implies £ = 0. Thus, is a phase function of positive 
type. Inspecting the partial derivatives of d^<f> w.r.t. x we conclude that the 
differentials d(d^ 1 4>), . . . , d(d^ n 4>) are linearly independent. □ 
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With a (z, .) € S°(X x R") we have exp[-Aa (z, .)] e 5°(X x R") by Lemma 
18.1.10 in 0. We define 

(2-9) 9(z',z)(x,0 ■= exp[-Aa (z,x,£)]. 

We shall keep this notation (for this symbol and others in the sequel) but it 
will be useful however to consider this symbol to depend on the parameters z 
and A instead of z and z 1 in the following analysis. Note that g( z >,z) is bounded 
w.r.t. z and c (o co w.r.t. A with values in S°(X x M"). Hence, we may define a 
distribution kernel G( Z ' ,z){ x> , x) € 9\X' x X) 

G(z',z){x',x) = J exp[i(x' - x\£,)]exp[-Aa(z,x',£)]$£ 

= J exp[i<j) { z',z){x',x,£,)} g( z ',z){x' ,0^ 

as an oscillatory integral. We denote the associated operator by G( z >, z )- This 
operator is often referred to as the thin-slab propagator (see e.g. jHlG])- We show 
that G{ Z ',z) is a global FIO in 1". 

Define a := (x',x, and 

u 9] (a,9) = d X] (j)(z<,z)(x',x,6) + £,• = -6j + £,-, 

u 6j (a,6) = d x ^(j) { z',z){x\x,e) = 6j -f< +iAd Xj ai(z,x',6), 

u X] {a,9) = d ej ^( z ', z )(x',x, 9) = x'j - xj + iAd^a^z, x', 9), 

where j = 1, . . . ,n. We denote by J( z ',z) the ideal in ^"^(R 5 ™) generated by 
the functions ue i , u^, and u Xj , and we let J( z ',z) be the subset of the functions 
in J( z ' ,z) that are independent of 9. 

Lemma 2.3. There exists Ai > 0, such that, for all z',z £ [0, Z], with z' > z 
and A = z' — z < A\ 7 the ideal J( Z ',z) * s generated by the functions 

(2.10) w (o) =d x ,.4> iz , tZ) (x',x,Z)-Z' j , 

= £j - £'j + iAd Xj a l (z, x', £) = u i} \ e = s , 
v Xj (a) = d£ j <j>( z i >g )(x',x,g) = x'j - Xj + iAd^ai(z,a/,^) = u X] \ e=i 

j = l,...,n. 

Some of the key arguments of the proof are close to that in the proof of Theorem 
25.4.4 in 0. 

Proof. The ideal J( z > jZ ) is also generated by the functions 

u 6j , Ufy := u 6] + = £j - ^ + iAd Xj ax(z, x', 9), u X] , 

j = 1, . . . ,n. We define v := (x',£',0), n := (x, £)• We set a point (^oMo) = 
(x' Q , £' , 9q,xo,!;o) where these generators vanish and we work in a neighborhood 
of this point. (Note that 0q = £o-) Since z i— > a\(z, .) 6 x K") is bounded 

we have that 3A X > such that for < A < Ai, and all z e [0, Z], 

det 9 , . . . , u e „ U£„ , u Xl u Xn ) jdv ^ 0. 
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By Theorem 7.5.7 in JOj we have 



(9)-(<r p( r , )(5) + (| ) ) 

where P is a < ^°° 2n x n matrix and Q is a < ^ , °° 2n x 2n matrix and the 
functions x, £ are also in a neighborhood of (uo,Ho). As the functions 
w x (u, n) := x' — x — x(/j,), w^(v, fi) := — £ — £(/i), wg{v, fi) := 6 — £ have linearly 
independent differentials, Lemma 7.5.8 in JU| proves that they generate J( z >,z) 
and the proof of that lemma shows that Q is invertible in a neighborhood of 
(vojfio). Letting 9 = £ we have 

ov,*,*,!*)- 1 ( W 'M ) - ( u ~f',' x i\ ) = ( v f\ ) . 

V w d v >N J V u d x > x >0 J \ v d a ) J 

We thus obtained that J( z ',z) is generated by the functions ug^ v Xj , v^, j — 
1, . . . , n. We then see that J( z ',z) is generated by v Xj , , j = 1, . . . , n. In fact, 
using Theorem 7.5.7 in ^U] again, any "if 00 function h(a) can be locally written 
in the form 

h{a)= ^2 (h X] (a',fi)v Xj (a',ti) + h i] (a',fi)v ij (a',fi)) + r(fi), 

\<i<n 

with a' = provided that < A < Ai. If h E J( Z '.z) then r E J( z >. z ) an d 

Lemma 7.5.10 in [TU] implies that VW E N, 3C N > 0: 

r(p) < Cjvmaxdlm x(jjl)\, |Im 

locally. Therefore, Theorem 7.5.12 in ^H] yields r E I(w x ,w^) = I(v x ,v^); 
which in turn implies g E I(v x , and thereby completes the proof. □ 

As the Poisson brackets (for the symplectic 2-form a' — a on T*(X' x X), where 
a' and a are the symplectic 2-forms on T*(X') and T*(X) respectively) of 
any two of the functions in (|2.1U|) vanish identically we obtain that the ideal 
generated by these functions is globally a conic canonical ideal in the sense of 
9, Definition 25.4.1. and Section 25.5]. The phase function 4>t z i tZ ) thus defines 
J( z ', z ) in the neighborhood of any point of J( z ',z)r'- it thus globally defines Jt z > lZ ), 
which is then of positive type. Therefore the operator G( z ',z) is a global FIO 
with complex phase (see Definitions 25.4.9. and 25.5.1. in 

Proposition 2.4. There exists Ai > such that if < A = z' — z < Ai then 
the operator G( z '. z ) is a global Fourier integral operator with complex phase and 
G (z , 2) eJ»(I'xl,(j (2 , 2) )',^ XI ). 

1 ll 

We denote the half density bundle on X' x X by ^x'xx- Note that ( J (z\z)Y 
stands for the twisted canonical ideal, i.e. a Lagrangian ideal (see Section 25.5 
in 0). 

Note that, with the following analysis, we could also consider Gi z ', z ) as a global 
FIO with real phase with amplitude in 5° (X' x X x JR n ) (see e.g. [201 )■ However 
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such consideration would be rather technical as one usually restricts oneself to 
the type S 1 ™ with p> \ for FIOs (see the remark at the end of Section 25.1 in 
see Also pages 391-392]). Viewing the thin-slab propagator G( z ',z) as a FIO 
with complex phase is also a good framework to understand the propagation of 
singularities in Part II. We shall however make this interpretation for G( z ',z) in 
Proposition 12 . 251 below, to apply a result of Kumano-go |131 Theorem 2.5]. 

We now establish some global continuity properties of the operator G( z ',z) stated 
in a slightly more general form (for similar results with global symbols see for 
instance |13| . where phase functions are real and other conditions are imposed 
on the phase function). 

Lemma 2.5. Let A be an FIO with a kernel of the form 

K A (x,y) = / exp[i<p.(x,0 - i(y\0]o- A (x,0<% £ @'(M. n x R"), 



where a A <E S m (R n xR n ) andip € ^°°(M n xM™) is such that Im(<p(x,£)) > and 
tp is homogeneous of degree 1 in for |£| large enough, and d Xi ip € S' 1 (R™ xR"). 
Furthermore, for all i = 1, ...,n we assume d^ i ip(x,£ l ) — Xi + fi(x,£) where 
fi G S°{W l x K"). Then A maps 5? into continuously. 



Proof. Let u £ 5? . We then have 

\Au(x)\< J \a A (x,0(i + \{;\)- m \\(i + \z\ru(0\<r(0 

<Csup \a A (x,d(l + \Z\)- m \ sup \(l + \^\r +n+1 u(0\, 

{eR™ £eR™ 

where C — J(l + |^|) _n_1 ^- The operator A is hence well defined from 5? into 
^(R™). If we differentiate we obtain 

D Xi Au(x) = J exp[ip(x,£)] (d Xi ip(x,£)a A (x,(;) - id Xi a A (x, £)) <f£ . 

Noting that d Xi <p(x,Z)a A (x,Z) - id Xi a A (x,Z) € 5 m+1 (M" x M") we similarly 
have 

\D Xi Au(x)\ < C sup |(1 + \t\) m+n+2 m\ 

CeR" 

< C sup \x a D p x u(x)\ for some a,/3 > 0. 
Iterating we find that Au € < ^ >co (W l ). Integrating by parts we also have 

A{xjU)(x) = J exp[i(p(x,0] (d^(p(x,^)a A (x,0 - id^a A (x,0) 
= XjAu(x)+ J exp[i<p(x,£)] (fi(x, i)a A (x, £) - id iz a A (x, 

Since /j^OM 35 * - ^M^, £ 5 m (R" x R n ) we obtain 

\x 3 Au{x)\ < C sup \x a D^u{x)\ +C sup la^'Df u(:e)|, 
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for some a, a',/3, /?' > 0. Similar estimates hold for \x a D@ Au(x)\ because of 
the hypothesis made on /j, i = 1, . . . , n. The operator A thus maps 5? into 5? 
continuously. □ 

To show continuity from 5^" into 5?" we shall need the following lemma. 
Lemma 2.6. Let j,k non-negative integers, u £ ,S fi {W rL ), f G C +1 (K") suc/i 

< Jm/(i) < C , iel", \f {r) {x)\ <C r , l<r<fe + l. 

TTien we have 



(2.11) ^' +fc 



u(x)(Im f(x)) J exp[iu)f(x)] dx 

<C^ SU P \D a u(x)\{\f'{x)\ 2 + Imf{x)) M/2 - k , uj > 0, 



|a|<fe xe 



where the constant C is bounded when the function f stays in a domain of 
'ta k+1 (W l ) where Cq, Ci, . . . , Ck+i can be chosen bounded. 

Proof. The proof is the same as that of Theorem 7.7.1 in ^U] where u G 1£q (M. n ). 
In fact the further assumptions on / made here allow to give global bounds that 
are needed since u G 5^ in the present case. □ 

Lemma 2.7. Let A be an FIO with a kernel of the form: 

K A (x,y) = /exp[z(x - y\0 + vyfat)] a A (x,0m G ^'(» n x R"), 



where a A G S m (R n x R") and 7 G S 1 (R n x R") is such that 1771(7(1, 0) > 0, 
and 7 is homogeneous of degree 1 in £, /or |£| Zarge enough. Furthermore, we 
assume that there exists d > smc/i t/ia£ 

(2.12) |i?e (9,7(2;, 0) I < d < 1, x G R", £ G R", |£| = 1. 

TTien A maps 5?" into ,5f" continuously. 

Observe that the differential of (f>(x, £) := (x — y\£) + 7(2:, £) does not vanish in 
R 2 ™ x R™\0. The function is thus a complex phase function. The differentials 
d(d^ 1 cf>) , . . . , d(d^ n 4>) are linearly independent. Hence is a non degenerate 
complex phase function of positive type. Note that by l|2.12|l the function (x — 
y\0 + li x i is an operator phase function in the sense of |5J Definition 1.4.4.]. 

Proof. Without loss of generality we may assume that 7 is homogeneous of 
degree 1 for |£| > 1. Let A f be the transpose of A and let u G S?, then. 

A t u{x)= J exp[-i(a;|f)] / exp[i(y|£) + ij(y, £)] a A (y, u(y) dySS, 

Define 

v(€,r))= / eKp[i(y\0 +ij(y,0}< T A(y,ri)u(y)dy , 
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and put w(£) = v(£,£). As u £ then v and w are both ff 00 . Then A*u is 
the Fourier transform of w. The lemma is proved if we show that u <— > w(£) is 
continuous from 5? to 5? . 

Let u = \i\ > 1 and Co - € We then have (y\0+l(y,0 = w/(l/,0 

with / homogeneous of degree in £, for |£| > 1. Note that d y f(y,f[) = £o + 
d y j(y,£o)- With the assumption made on d y ^y we have \d y f(y,£)\ > c > 0. 
Applying Lemma f2. 61 and estimate (|2.11() we obtain 

\a\<k 

<K' k (l + \ V \) m sup \D a u(y)\, u>l 

\a[<* 

where the constants Kk, K' k can be chosen uniformly w.r.t. £, |£| > 1 since the 
constants Co, Ci, . . . , Cfc+i of Lemma l2. Cl ean be chosen bounded (as £o € 
Now setting 77 = £ we obtain that for all fc € N, 3ift > 

(2.13) (1 + |e|) fe -"XOI < ^fe sup |£><XlOI ,Cer, |£| > 1. 

I a\<h 

We now consider 

D &H0 = / ex p[«(ylO + *7(#>f)] 

((y* + %7(2/,0) cr A(y,0 - id ii a A (y,0) u (y)dy ■ 

As yiu{y) G =5^ and %7(y,£) is homogeneous of degree for |£| > 1 estimates 
similar to those in (|2.13(l are valid. □ 

It is immediate from the structure of (f>( Z ',z) m 1|2.8[) that Lemma 12.51 applies to 
S(z',z)- If A = z' — z is small enough we have \Ad Xi bi(z,x',£ l )\ < d < 1, due to 
Assumption ll.il and then Lemma \l . 71 applies . We thus have 

Proposition 2.8. There exists A2 > such that if z\ z € [0, Z] with < A := 
z' — z < A2 then G( z ',z) maps & into ,5^ and S^" into ,5f" continuously. 

Remark 2.9. By the above result, composition of the two FIOs G( z " tZ ') an d 
S(z',z) is thus natural without further requirement such as having the operators 
properly supported. 

We now turn to global L 2 and Sobolev space continuity for the operator G( z ',z)- 
We shall use the following lemma. 

Lemma 2.10. Let p s (x,t;) be bounded w.r.t. the parameter s with values in 
S™(X x R") and define 

£(A,x,€) :=Z-Af(x,0 

where f is in S 1 (X x R",R") and homogeneous of degree 1 in £, for |£| > 1. 
Then 

Ps(A,x,0 :=p s (x,£(A,x,£)) 
is bounded w.r.t. s with values in < ^' oo ([0, A3], S™(X x R™)) for A3 small enough. 
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Proof. Let A 3 be small enough such that |£ — A/(x, £)| > Co > if |£| = 1 and 
A € [0, A 3 ]. We then have 

i + Coiei < l + IC-A/^OI <i + Ci|e|, £eR", |£| > 1, Ae [0,A 3 ]. 

This inequality yields the proper estimates for d"d^p s to prove that p s € 
S™(X x R n ). As derivatives w.r.t. A do not affect the symbol order and type, 
the proof is finished. Bounds w.r.t. the parameter s come naturally. The proof 
is complete. □ 

Following [22 we introduce 

Definition 2.11. Let L > 2. A symbol q(z, .) bounded w.r.t. z with values in 
S 1 (MP xl r ) is said to satisfy Property \Pl\) if it is non-negative and satisfies 

(P L ) \^q(z,y, V )\ < C(l + \r,\)-W+(\<*\+WV L 

(l + q(z,y,r } )) 1 - {laMm/L , ze[0,Z], y eW p , v eR r . 

We then set p = 1 — 1/L and 5 = 1/L. 

Remark 2.12. Suppose q(z, .) as in Definition 12 . 1 II and \a\ + \(3\> L then 
(1 + | ?7 |) 1 -(I«I+^I)/- L < C(l + q(z, y, jjJjMM+M)/* z G [0, Z] , y S R p , 77 e R r . 
Estimate jP^ is thus clear in this case. 

Examples of symbols with such a property with L > 2 are given in |21| . In fact 
we prove that c\ satisfies Property \Pl\ for L = 2. 

Lemma 2.13. Let q(z,y,r/) be bounded w.r.t. z with values in S 1 (M. P x R r ). If 
q > then q satisfies Property \Pl\) for L = 2. 

Proof. Bounds w.r.t. z are natural; we shall omit the dependence on z in the 
proof for concision. We have to prove that 

\d«dPq\ < C (1 + |r7|)5(M-l/3|) (1 + 

The property is clearly true for \a \ + \f3\ = and for \a\ + \ (3\ > 2 by the remark 
above. Let us now treat the case \a\ + \f3\ = 1. For this we use Landau's 
inequality: let / G C 2 (R) with / > and /" is bounded then (see [U page 40] 
and HOI Lemma 7.7.2]) 

|/'(*)|< 2 (/(*))* fsupirwi 

\teR 

We first treat the case \a\ = 1. Define p(y,r)) = (1 + M 2 )~^ q(y,r/). Then 
p € x R r ) and <9 2t >(y, 77) is in S* (R P x R r ) and is thus bounded. We thus 

have 

(1 + M 2 )^ l^(y,ry)|<C((l + M 2 )^ q(y, V ))K 
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which yields 

\d^q(y,ri)\<C(l + \ri\)i (1 + g(|/, »/))*, 

which is the expected estimate. Let us now treat the case = 1, with for 
instance, /3 = (1,0, ... ,0) and a = (0, . . . ,0). Define p(y, rf) = (1 + |?/| 2 ) 5 q(y,v)- 
Then p e S 2 (W x R r ) and thus d^p{y, rf) is bounded. We hence have 

\dfo(y,r))\ < C(p(y,r)))i. 

With 

d^p(y, V ) = (l + \ V \ 2 ^ d^q(y, V )+ m (l + \r ] \ 2 )-^ q(y,r,), 
the triangular inequality yields 

(i + M 2 )* \tfq(y,v)\ <C( P (y^ + \vi\(i + H 2 r h q(y,v) 

<C(q(y, V ))^(l + \nf)i + (q(y,n))^ 

<C{q{y,r,))* (1 + |„| 2 )*. 

We finally obtain 

\tfq(y,v)\ <C(q(y, V ))H^ + \v\rK 
which is the expected estimate. □ 

Remark 2.14. If the symbol q(z,y,rj) satisfies Property ]Pl$ then the ampli- 
tude q(z, y 1 , 77) + q(z, y, rf) also satisfies Property \Pl\ (with derivatives w.r.t. y, 
y' and rf). 

Proposition 2.15. Let q(z, .) be bounded w.r.t. z with values in S 1 (MP x M r ) 
with q(z,.) > 0. Let q(z,.) satisfy Property \Pt\) and define p/\{z : y,rf) — 
exp[— Aq(z, y,T})]. Letm^z N. Thenq m p/\ is smooth w.r.t. A, bounded w.r.t. z, 
with values in Sp(R p x M r ) for A in any interval [A m j„, A max ] with A m i n > 0. 

Proof. dydP(q m pA) is a linear combination of terms of the form 

A k (d^d b Jq) . . . (dpd h v 'q)(d^d^q) . . . (d^d^q™- 1 p A 

with < I < m and a%+- ■ -+ai+ai + . . . = a and b\ + - ■ •+6j+/3i + . . . (3k = ft- 
We can estimate the absolute value of each of these terms, using Property \Pl\ - 

by 

CA fc (l + 1771)^1+^ (l + g)<+ fe -^g™-< M 

as (1 + (? )'+ fc -^T^g m - z A fe + m / 9 A < C. □ 

While the symbol exp[— Aq(z, y, rf)] is bounded w.r.t. z and smooth w.r.t. A 
with A > A min > with values in S°(W x W), this fails to be true at A = 0: 

9 A ex P [-Ag]| A=0 = -qi S° p (W x E r ). 

In fact when we want to control the behavior of exp[— Aq] close to A = we 
shall use the following definition and lemmas. 
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Definition 2.16. Let L > 2, p = 1 — 1/L and 5 = 1/L. Let pA{z,y,n) be a 
function in %f°°(R p x W) depending on the parameters A > and z G [0, Z]. 
We say that p<\ satisfies Property \Ql\) if the following holds 

(Ql) aX(PA- PA\A=o)(z,y,v) = A m+5 W+Wp™ a(3 (z,V,r)), 

for \a\ + \{3\ < L, < m < 1 - S(\a\ + \/3\), 

where p T ^ a,3 (z,y,ri) is bounded w.r.t. A and z with values in 5™ p l^l +5 l a ! (RP x 
R r ). Lt follows that pA(z,y,n) — pA\A=o{z,y,r]) is itself bounded w.r.t. A and 
z with values in S°(W xW). 

Lemma 2.17. Let q(z,.) be bounded w.r.t. z with values in S 1 (MP x M r ) and 
satisfy Property \Pl\) . Define pA(z,y,r/) = exp[— Aq(z, y, rf)]. Then pa satisfies 
Property \Q L [ for A S [0, A max ] for any A max > 0. As pa|a=o = 1, PA is 
itself bounded w.r.t. A and z with values in Sp(W p x R r ). 

Proof. In the proof all the functions and symbols will naturally be bounded 
w.r.t. z. We thus drop the variable z here for concision. 

We define 

PT P ■■= A- m - 5 (H+l/3|)^ (jOA _ PA | A=0 ). 

We first consider the case \a\ + \(3\ — with < m < 1. We need to estimate 
\dyd b p A l0a \. The case m = 0, \a + b\ = has to be treated independently 
but is trivial: we clearly have = \pa — 1| < C. We shall now estimate 

I<9^Pa°°I = \A- m d$d b (p A - 1)| in the case where m > or \a\ + \b\ > 0. For 
this we write 

(2.14) p A (y,r)) -1 = -A q(y,T])exp[-sAq(y,r))]ds. 

Jo 

We then have Pa°°(2A V) = — Jo 1 Qa( s i Vi r ?)^ s with 

1a(s,V,v) = & 1 - m q(y,r])exp[-sAq(y,ri)]. 

We prove that 

\d a y d b v qt{s,y,ri)\ < C(s)(l + |,,|)*»-/>l*l+*l«l 
with C(s) bounded w.r.t. A and L 1 w.r.t. s e [0, 1]. The result then follows for 

„m00 

Pa ■ 

When computing dyd b q A we obtain a linear combination of terms of the form 

A l - m {d a y a d^q){-sA) k {d^d h r /q) . . . (d^d^q) ex V [-sAq], 

with a + a 1 + ■ ■ ■ + a k = a, b + bi + ■ ■ ■ + b k = 6, 

(where k can be 0) . Using Property \Pl\ we find that the absolute value of such 
a term is bounded by 

CA 1 - m (sA) k (l + \ r} \)-\ b \+ s <.\ a + b \)(l + q) k+1 - S ^ a+b ^ e^[-sAq] 
< Cs m+5 ^ a+b ^^ 1 (l + |^|) m_ ' , l b l+' 5 l a l A 5 ^ a+b ^ 

(sA(l + q ))-m+k+l-S(\a+b\) exp [_ sA? ] j 
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as 1 < C(l + M) m (l + q)~ m if m > 0. If I := -m + k + 1- S(\a + b\) > wc 
use that (sA(l + gr)) z exp[-sAg] < C if < s < 1, < A < A max and q > 
and we obtain the following estimate 

^ s m+i5(|o+6|)-lQ ^_ |^|yn-p|h|+5|a|^i5(|a+6|) 

If Z < 0, (1+q) is simply bounded (q > 0) and we obtain the following estimate: 

CA k+1 ~ m s k (l + | 7 ;|) m -' > l fc l+ ,5 l a l. 

As m + 5(\a + b\) — 1 > — 1 in the considered case, both estimates exhibit bounds 
that are in i 1 ([0, 1]) w.r.t. s. We also have uniform bounds w.r.t. A as we have 
assumed m < 1. 

We now treat the case 1 < \a\ + |/3| < L, < m < 1 — 5(\a\ + \/3\). We estimate 
the absolute value of 

dyd^pT ) = A- m - 5( ^+W d« +a d b + f3 PA 
which is a linear combination of terms of the form 

A k-m-6(\a\+m (d^d^q) . . . (d%» d b J q) exp[- Aq] , 

with <2i + • ■ ■ + at = a + a, b\ + • • ■ + bk = b + 0, 

where k > 1. Using Property ]PlI we find that the absolute value of such a 
term is bounded by 

CA k-m-S{\a\+\0\)r 1 + t v \\-\0\-\b\+6(\a\+\a\+\0\+\b\) 

(l + ? )*-*<M+l»l+l0l+Wexp[-A«] 
< C(l + | 7? |) m -p(l' 3 l+l fc l)+ ,5 (l Q! l+l a l)(i + q y s (\ a \+\ b \) 

(A(l + q ))k-m-8(\ a \+\0\) exp [_A g ] 

< C(l + |^|)™ l_ ' , (l' 3 l+l fc l)+ ,5 (l Q l+l a l) 5 

as k - m - 5(\a\ + >l-m- 5{\a\ + \(3\) > and < A < A max . This 
completes the proof. □ 

Lemma 2.18. Let f E Sf°°(M) and q A (z,y,r]) m <g°°(MP x M r ) that satisfies 
Property \Ql\1 arid such that <7a(z, .)|a=o is independent of y and r\. Then 
f(q A )(z,y,r)) satisfies Property \Q L \> - 

Proof. Again bounds w.r.t. z are clear. We first treat the case \a\ + \(3\ = 0. We 
write 

/(<Za) - /(<Za|a=o) = (?A - <7a|a=o) / /'((I - s)(7a|a=o + sq A )ds. 

Jo 

As <?a|a=o is independent of y and rj, then q A is bounded w.r.t. A with values 
in S°(W x W) by Property j^IJ and so are (1 - s)q A \ A=0 + sq A and /'((l - 
■s)<7a|a=o + s 9a) by Lemma 18.1.10 in [5] with bounds in 5"(R P x K r ) uniform 
with respect to s. We thus obtain that J„ /'((l — s)9a|a=o + sq A )ds is bounded 
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w.r.t. A with values in S°(R P x W). We conclude using Property \Ql\ for 
9A — <Za|a=o- Let us now treat the case 1 < \a\ + \f3\ < L and choose < m < 
1 — S(\a\ + We see that dyd@f(q A ) is a linear combination of terms of the 
form 

(d^d^q A ) . . . (d^d^q A )f^(q A ), 

where k > 1, a\ + ■ ■ ■ + = a, j3\ + ■ ■ ■ + j3k — (3. Now choose < < 
1 — <5(|cti| + i — 1, • • • , k, such that m — mi + • • • + rrifc. Then Property 
J[\ yields terms of the form 



^TOi+5(|ai| + |/3i|) fcm k +8{\a k \ + \P k \) qmiaiPi ^m k a k k _ ^m+<5(|a| + |/3|)^ma/3 

with q r ^' a% ^\ i = 1, . . . , k, bounded w.r.t. A with values in s™'~ p ^ a ^ +S ^^ (W x 
W) and q^" := q^ 1 "^ 1 ■ ■ ■ q™ kak/3k . We note that /« {q A ) is bounded w.r.t. A 
with values in S^R^ x R r ). The symbol g^ Q/3 is bounded w.r.t. A with values 
in s 7 p n - plal+sm (W x R r ), which yields the result. □ 

With Remark OH Lemma EH£1 and the previous lemma we obtain 

Corollary 2.19. Let f 6 < if 00 (R) and let q(z, .) bounded w.r.t. z with values in 
S 1 (W x R r ) satisfy Property JF^p. Define 

PA{z,y',y,v) = exp[- A(q(z, y', rj) +q(z,y,r)))]. 

Then f(p A ) satisfies Property \Ql\)- As /(j»a)|a=o = 1 ; I{pa) is itself bounded 



w.r.t. A and z with values in 5"(M 2p x R r ). 



Note that the property \Ql\ is stable when we go from amplitudes to symbols: 
Proposition 2.20. Let q A (z,x,y,£) be an amplitude !nSj(l 2p xl p ) depending 



on the parameters A > and z € [0, Z] that satisfies Property Then 
G {qa} (2, x. £) satisfies property \Ql)I- 



Proof. We use the oscillatory integral representation for the symbol: 

a {q A } (z, x, £) := / / exp[-i{y\r))} q A {z, x,x - y,£ - rf) drj dy. 



Let < |a| + \/3\ < L and < m < 1 - S(\a\ + |/3|). Computing d^d^{a{q A } ~ 
g {qa} |a=o)j we obtain a linear combination of terms of the form, with a±+a2 = 
a. 

J J exp[— i{y\r)}] d^d^ 2 df (q A - q A \ A =o){z, x, x - y, £- 77) otr] dy 

= JJ exp[-i(y\r,)} A m + 5(\ a \ + \0\) q M«u**W { ^ x ^ x _ y ^_ v)Srt dy 

_ ^m+S(\a\ + \p\) a . ^m{ ai ,a 2 )f3^ 

where q™^ 1 ' 012 ^ [ s bounded w.r.t. A and z with values in the symbol space 
Sp l p l /3 l +(5 l Q! l (K 2 p x R p ). As the map a <— > a {a} maps bounded sets into bounded 
sets the result follows. □ 
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We shall also need the following lemma. 

Lemma 2.21. Let q A (z,x,y,£) be an amplitude in 5°(M 2p x R p ) depending 
on the parameters A > and z G [0, Z] that satisfies Property \Ql\1 for 1 < 
M + < 2 and such that q A (z, .)|a=o * s independent of (x,y,£). Let r(x,£) G 
S S (W x RP) for some s G M. Then 

a {q A r} (z, x, - q A (z, x, x, £) r(x, = A m+2S \%(z, x, £), < m < p - 5, 

where the function X l A (z,x,^) is bounded with respect to A and z w»t/i values in 
S™ +s - {p - S) {W xRP). 

Proof. For the sake of concision we take p — 1 in the proof but it naturally 
extends to p > 1. We write Aa = q A r. Using the oscillatory integral represen- 
tation of a {q A } we obtain 

<r{q A r} (z,x,£) - q A (z,x,x,£)r(x,£,) 

Taylor's formula yields 

o-{q A r} (z,x,£,) - q A (z,x,x,£,)r(x,£,) 

-yexp[— i{y\£-r])] d 3 X A (z,x,x- sy,rj) tin dy ds. 
With an integration by parts we obtain 
o-{q A r} (z,x,£) - q A {z,x,x,£,)r(x,£,) 



i) 



iexp[— - n}] d$d±\ A {z, x,x— sy,n) &r\ dy ds 

= a\i dsdi\ A (z,x, (1 - s)x + sy,£) ds 



i) 



I Jo 

where d 3 d 4 X A (z,x,y,^) = {d y d^q A ){z,x,y^) r(x,£) + d y q A (z,x,y,£,)d i r(x,£,), 
as r does not depend on y. The first term is treated using Property \Ql\ while 
for the second one we write 

a a \m' +8 m' (0,1)0 o 

O y q A <%r = A m + q A d s r, 

where < ml < 1 - 5 and oX' ( °' 1)0 r G S 7 / +s - 1+s (M. 2p x M p ) by Property tfJJ}. 
We actually take £ < ml < 1 — 6 and write m = ml — 8. We obtain 

d v q A d^r = A m+2S ft, 

where q% is bounded w.r.t. A with values in S™ +s -p+ 5 (M?p x W) and < m < 
1 — 2(5 = p — 5. We conclude since the map cr{.} maps bounded sets into bounded 
sets. □ 

We are now ready to give an estimate of the operator norm of the thin-slab 
propagator, </(*>,*), in L(H^ {X), {X')) for any set. 
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Theorem 2.22. Let set. There exists M > 0, A 4 > such that 

\\Sv,z)\\(HM,HW) < 1 + AM, 
for all z', z G [0, Z] s-uc/i that 0<A = z'-z<A 4 . 

In the proof we assume that C\ satisfies property Pl for some L > 2. We 
know that it is always true for L — 2 by Lemma 12.131 but special choices for 
Ci can be made. As before we use p = 1 — l/L and S — 1/L with p > 5 
for L > 2 and p = <5 = ~ for L = 2. In the proof we proceed classically by 
computing G( z ',z)Sf z ' z ) an d use the classical results on -0DOs (see e.g. ^0 
Section 5] and also 0). Here we however do not content ourself with the 
continuity of G( z >, z ) but we want to obtain a precise estimate of the operator 
norm in L{H^ (X), (X')), which will be required in the sequel. Here we 
exploit the fact that A can be taken arbitrarily small which allows to carry out 
some explicit computations. 

Proof. Let s G M, then the kernel of Ar z ' :Z ) := G( z ', z ) ° E^'^ is given by 

•A( z > tZ )(x',x) = J exp[i<j>( z , iZ) (x',x,£)] ff( 2 ', z )(x',£) (£)~ s <f£. 
Computing the kernel -D(v. z ) of 2?( z ', z ) := A( z > iZ ) ° A* z , x * we obtain 

= J exp[i(x' - x\£) +iA(bi(z,x',g) - h(z,x,£))]d( z i tZ )(x',x,g) <f£ 

where 

d( z >, z )(x' ,x,C) 

= exp[-A(c\{z,x' ,£) + c x {z, x, (x' , £,) g( z >, z )(x,£) (0~ 2s - 

We write b\{z, x', £) — bi(z, x, £) = (x' — x\h(z, x', x,£)) where h is smooth and 
homogeneous of degree one in £, |£| > 1. The function h and continuous w.r.t. z 
with values in S 1 (X x E n ) by Assumption II. II and estimate (1.1.9) in [Ej. We 
thus obtain that the change of variables £ — ► £ + Ah(z, x', x, £) = -ff(A,z,z',a;)(0 
is a global diffeomorphism for A small enough (uniformly in z G [0, Z]). We 
denote £(A, z, x', a;, £) = H7 Azx , x j(0- We thus have 

D (z',z){x',x) = J exp[i{x' - x\£)]d( z > tZ )(x',x,€(A, z,x',x,£)) Ja(z,x',x,£) <f£ 
where Ja(z, x', x, £) is the Jacobian. 

Lemma 2.23. T/ie function £(A,z,x',x,£) is homogeneous of degree 1 in £, 
/or |£| > 1, continuous w.r.t. z, H 00 w.r.t. A wifft values in S 1 (R 2n x (R n )) if 
A is small enough, i.e., 

3A 4 > 0, £e tf°([0,Z],^([0,A 4 ],S 1 (K an x (IT)))). 
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Proof. Homogeneity is clear. We have 

\i(A,z,x',x,£)\ = \£-Ah{z,x',x,£(A,z,x',x,£))\ 

<l + AC(l + \i(A,z,x',x,0\), |£| = 1, 
which yields, because of homogeneity, 

\z(a,z,x\x,o\<±±^(i + \z\), iei>i, 

for A small enough, uniformly chosen w.r.t. z £ [Q,Z], x',x £ R n . Differentiat- 
ing the j th coordinate of £, 

£? = £7(^,2, + Ahj{z,x' ,x,l{A,z,x' ,x,C)), j = 1, ■ ■ ■ ,n, 
w.r.t. Xi yields 

(2.15) d Xi £j(A, z, x', x, £) + Ad^hjiz, x', x, |(A, z, x', x, £)) 

+ A ^2 d i l h j( z ^ x '' x ^(^,z,x',x,^)) d x Ji(A,z,x',x,€) = 0, 
1 

j = 1, . . . ,71. 

The partial derivatives of h are bounded for |£| = 1. We can hence solve for 
3 Xi £(A, z, x' , x, £) when A is small enough and find the expected estimate from 
that obtained for £(A, z, x' , x, £): 

3O0, \d Xi i(A,z,x\x,0\ <C(1 + |£|), a^areR", ( e 1™. 

Differentiating w.r.t. x[, £j, and A yields similar structures and the proper 
symbol estimates. The proof carries on by induction. Note that the required size 
for A to solve the systems of the form (12.15(1 remains fixed along the induction 
process. □ 



Continuation of the proof of Theorem \2.22\ From (the proof of) Lemma l2\23l we 
also obtain that the Jacobian Ja{z, x\ x, £) is homogeneous of degree zero in £, 
|£| > 1, and is continuous w.r.t. z and < ^' 00 w.r.t. A with values in S^R 2 ™ x W 1 )). 

We write Pa{z, x 1 , x, £) := exp[— A(ci(z, x',£) + ci(z, x, £)]. As ci satisfies Prop- 
erty \Pl\ we then have pa satisfying property \Ql\ by Corollary 12. 191 Define 
Pa(z, x', x, £) := pa(z,x',x,£(A,z,x',x,£)). Lemma fe.lOl and Lemma r2.23l vield 
that pa satisfy property \Ql\ ■ We then have 

d( z >,z)(x\x,£(A,z,x',x,£)) J A (z,x',x,£) =:p A (z,x',x,£,) k A (z, x' , x, £) 

where Ua(z, .) is bounded w.r.t. z and c /a°° w.r.t. A with values in S~ 2s (M. 2n xM. n ) 
and k A {z, .)|a=o = (-)~ 2s by Lemma l2~Tni and Lemma l2~2"31 By Theorem 1.1.9 
and formula (1.1.9) in JHj we obtain 

k A (z,x',x,£) = (0~ 2s + Ak A (z,x',x,0, 

where Ik a is bounded w.r.t. z and < ^'° w.r.t. A with values in S~ 2s (R 2n x l n ). 
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Call T [z ^ z) := £W o V {z t, z) o Its symbol is in S° p {W l x R") and is given by 

/(*,«)(*', == ((0 s # <T{p A (2,a/,a:,0 M^s.O} # (O'Xs'.O 
= ((0 s # a{p A (z,x',x,Q(0- 2s } # (0 S )(^',0 

+ A((0 S #a{p A (z,x', a; ,0fcA(^x'^,0(0" 2s } # (O*)(*',0 

As pa bounded w.r.t. z and A, A small enough, with values in Sp(M. 2n x K n ) 
(Property \Ql\ ) we obtain that the second term in the equation above satisfies 
the same property and thus we can write 

where T^ z , ^ has for symbol 

((0 S #*{PA(Z,X',X,Z)(1;)- 2S } # (0 S )(Z',0 

and \\^ z ' Z ^\\(L 2 ,L 2 ) 5~ -K" ) uniformly in z G [0, Z] and A, A small enough, 
by the Calderon-Vaillancourt theorem (see ^] Chapter 7, Sections 1,2] or jJSJ 
Section XIII-2]) in the case L = 2 and by Theorem 18.1.11 in [H] in the case 
L > 2. With Lemma 12.211 we see that 

a {pa (*, x', x, (0 - 2s }-pa(z,x',x',0 (£) ~ 2s = AA A (z,x\ 

where Aa is bounded w.r.t. A and z with values in S~ 2s (M. n x K n ). We thus 
obtain 

•F(z',z) = ^{z',z) + ^F{z',z) 
where T h ^ z , z \ has for symbol 

/i(z,x',0 :=({0 S #Pa(z,x',x',Z){0- 2s # (0 S )(^',0 

= ((0 S #PA(z,x',x',0(^r s )(z,x',0 

and \\J-"? Z , z \\\(l 2 .l 2 ) < uniformly in z G [0, Z] and A, A small enough. 

For the rest of the proof, if we don't write it explicitly, by pa and pa{z, x,£) we 
shall actually mean Pa(z, x,x,£). 

Lemma 2.24. 

((•) s # PA{z 7 .)(.y s )(z,x,Z) - PA (z,x,£) = A(j, A (z,x,£), 
where [Xa(z,x,£) is bounded w.r.t. z and A with values in S®(X x R n ). 

Proof. We write 

Pa(z,x,£) = (0 _s // exp[— i{y\£ - rj)] (r]) s pa{z,x,£) dr] dy 
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and thus obtain, with the oscillatory integral representation for the composition 
formula, 



(0~ s // exp[-i(y|£- r))} (r)) s (p A (z,x - y,£) - p A (z,x,£)) d~r) dy. 



With Taylor's formula and applying an integration by part, we find (we have 
supposed n = 1 for the sake of simplicity but it naturally extends to p > 1) 

((.) s # Pa (z,.)(.)- s )(z,x,0-Pa(z,x,0 = 

(0 _s J J J «cxpH(y|£ - 77)] d,-!^)* d x p A (z,x- ry,£) ctrj dy dr. 

Using Property \Ql\ with m = 1 — 5 we find 

((.y # PA(z,.)(.)-)iz,x,t)- PA (z,x,Q = A(t)- S 

iexp[-i{y\Z-r))] d,^!]) 8 q A w (z, (1 - r)x + r(x - y), f) <fr? dy dr 

= A(0- s (d ( (O s # q^ 10 (z,u,x,0)\u^ 

where 

q% w (z,u,x,0= [ q% w (z,(l-r)u + rx,0 dr. 
Jo 

As q 7 ^ 10 is bounded w.r.t. A and z with values in S*(R ra x R n ) we obtain the 
result. □ 

End of the proof of Theorem \2.22\ With the previous lemma we see that 

where !F? Z , g s has for symbol p A (z, x, x' , £) and \T^ X , z ^\\(l 2 .l 2 ) < K 3 uniformly 
in z G [0, Z] and A, A small enough. 

To estimate \\J-^ Z , ^ || (l 2 ,l 2 ) we follow the procedure at the end of the proof of 
Theorem 18.1.11 in 8 1. Let A := 1 + A. Define 



v A {z,x',£) = V 'A 2 - \ PA {z,x',x',£)\ 2 , 
which satisfies Property \Ql\ by Lemma \2. 181 Then define r A by 
v A jfv* A = A 2 -p A jfp\ - r A . 

Note that 

v A # a?', f) = a {v A (z, x', £)vK{z, x, £,)} (z, x, £) 



It is easy to check that v A (z,x', £)v A (z, x, £) satisfies Property JQlJ for |a| + 
> 1. The same applies to p A (z,x',S,)p A (z,x,^). Lemma 12.211 applies and 
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with m = p — 8 we thus obtain that ta = Af a with f/\ bounded w.r.t. z and A 
with values S° p {X xl"). Thus 

||^ V )||(^) = < + A) 2 + AC < 1 + AK 4 , 

for some K 4 > large enough. We thus obtain that \\J-( Z ' , z )\\(l 2 ,l 2 ) < 1 + AK 
where K = K 1 + K 2 + K 3 + K 4 . With the definition of T {z ^ z) it follows that 

\\G(z',z)\\(m°l,H<.°i) = \\(G(z>,z))*\\(hw,hW) ^ VTTAK 
which concludes the proof of Theorem 12. 221 □ 

We observe that for A small enough, the function (x'\£) + Abi(z,x',£) satisfies 
the conditions (P)-(i), (P)-(ii), and (P)-(iii) in [I'M page 2]. With Lemmas l2.13l 
and 12.171 we observe that an FIO with phase function 4>( Z ',z)( x '> x i an d am- 
plitude in <ja(z, x' , £) in S m (X x R) may actually be understood as an FIO with 
real phase (x' — x\£) + Abi(z, x',^) and amplitude (Ta(z, x', £) exp[— Ac%(z, x 1 , £)] 
in S™{X x R). Applying Theorem 2.5 and the following remark in we obtain 

Proposition 2.25. Let -4( z /. z ) be the global FIO with kernel 
A(z',z){x ,x) = J exp[i<f>( x , i!s )(x',x,g)] a a(z,x' 

with o~a(z, .) bounded w.r.t. z with values in S m (X x R n ), m £ R. TTien /or all 
set t/iere exists M = M(s, m) > 0, A 5 > smc/i iftat 

\\^(g',z)\\(HW,m-^)) < M p(va(z, .)) 

for all z £ [0, Z], and < A < A$, where p(.) is some appropriately chosen 
semi-norm in S m (X x R"). 

This proposition could also be proved by adapting the proof of Theorem 12.221 
to this case. Note that in the case a a = 9(z',z) we were able, in the proof 
of Theorem 12.221 to achieve a finer estimate. The proof heavily relies on the 
particular structure of the phase function and the amplitude that can be taken as 
"close" as we want to those of the identity operator by taking A small enough. 



3 The approximation Ansatz. Convergence in 
Sobolev spaces 

We first define the Ansatz that approximates the solution operator to (|1.5|l - 
(|1.6(l . We chose to use a constant-step subdivision of the interval [0, Z] but the 
method and results presented here can be naturally adapted to any subdivision 
of [0,Z]. 

Definition 3.1. Let ?p = {z^°\ z^ x \ . . . , z^} be a subdivision of [0,Z] with 
= z(°) < < • • • < z<W = Z such that zS i+ ^ - z& = A<p. The operator 
Wq3 jZ is defined as 

{ G(z,o) if Q<z<zW, 



i=i 
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The following proposition will be useful in the sequel. 

Proposition 3.2. Let s G K. There exists K > such that for every subdivision 
<P = {z^, . . . , z^} of[0,Z] with = z(°) < < ■ ■ ■ < zW = Z and 
Wtp jZ as defined in Definition \3.1\ we have 

We[0,£], ||W !p ,,|| (ffW ,ffW)<-K', 
if Asp is small enough. 

Proof. By Theorem 12.221 there exits M > such that if A = z' — z is small 
enough then ||{?( z ', z )||(ijoo jjwj < 1 + AM for all z € [0,Z]; we then obtain 
l|VW«p, a ||(u-(.),fr(«)) < (l + AqjA/)^ = (l + ^r)^ which is bounded as it converges 
to exp[ZM] as N goes to oo. □ 

It should be first noticed that Wyt :Z is not the solution to problem H1.5() - (|1.6() 
even in the case where the symbols b and c depend only on the transversal vari- 
able, (x). While singularities propagates along the bicharacteristics associated 
with ia\ — 6i, observe however that, with the form of the phase function 4>( Z ' ,z) 
in 12.8fl . the operator G( Z ',z) propagates singularities along straight lines. See 
Part II, for further details, in particular regarding the set Jr z ',z)R that replaces 
the canonical relation for the propagation of singularities for FIOs with complex 
phase |H1 Sections 25.4-5]. 

Furthermore, by composing the operators Gt z ",z') an d Q( Z ',z)i °ne convinces 
oneself that 

G(z",z) 7^ G(z",z>) ° G(z',z) 

in general if z" > z 1 > z G [0, Z] (use again that singularities propagate along 
straight lines). The family of operators {G{ Z ' ,z)) (z 1 ,z)e[o,Z] is thus neither a semi- 
group nor an evolution system. 

We now proceed towards the proof of the convergence of W<p jZ to the solution 
operator to problem l|1.5|l - i|1.6[l m the sense of Sobolev norms as N = |$p| goes 
to oo. 

Lemma 3.3. Let set and z",z G [0, Z], with z < z" . The map z' i— > G( z > lZ )> 
for z' G [z, z"\, is Lipschitz continuous with values in L{H^ S+1 \X) 1 H^ S '(X)), 
for z" — z — A small enough. More precisely there exists C > such that for 
all u Q G H( S+1 \X) and z«, z^ G [z, z"] 

(3.16) || - G izm , z) )(u )\\ Hls) < C\zW - z^\\\u Q \\ H(s+1) . 

Proof. Let z^\z^ G [z",z] and let u G H S+1 (X). Write 

(G(zW,z) - G {z {i)^)){uo){x') = 

/ / exp[i(x' — x\£) — (z — z)a{z, x , £)] a{z, x , £) uq(x) dx ct£ dz . 
(i) J J 

When A is small enough we can apply Proposition 12 . 251 and obtain H3.16(l □ 
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Lemma 3.4. Let s € R, z",z G [0,Z], wit/i z < 2", and Zet m € # (s+1) (X). 
TTien ifte mapz' i-> ^, 2 )(tt ) «s m tf°([z, z"}, H^ +1 \X))f\^ x {[z , z"], ffW(l)) 
/or z" — z = A small enough. 

Proof. Let G [z, z"] and let e > 0. Choose z" — z small enough such 
that Theorem 12.221 and Lemma 13.31 apply and Choose u\ G H ( s + 2 ) such that 
\\uq — Ui\\ H (s+i) <e. Then for z^ G [z,z"] 

(3.17) \\G {z (2) tZ) (u ) - G {z (i) tZ) (u )\\ His +i) < \\G( z (2) >z) (u - 

+ \\G(zV) ,z)( u l) - G( z (i) tZ )(ui)\\ H ( s + i) + \\G( z (l) <z )(u - Ui)|| ff (,+i) 

< 2(1 + AM)e + C\z^ -« (1) |||«i|| H (.+ a) . 

The continuity of the map follows. Differentiating Gr z t tZ \(uo) w.r.t. z' we can 
prove that the resulting map z' 1— > d z >Gt z i <z ){uq) is Lipschitz continuous with 
values in L(H^ s+2 \ H^) following the proof of Lemma f3. 31 there exists C > 
such that for all v G H^+ 2 ^ (X) 

\\{d z ,G {z v), z) - d z ,G {z a), z) ){v)\\ HM < C\zW - zW\\\v\\ Hl . +a) . 

We also see that the map v 1— ► d Z 'Gi z i jZ \{v) is continuous from Jf( s+1 ) into 
H ^ with bounded continuity module according to Proposition 12.251 With 
uq G i/' s+1 )(X) we make a similar choice for Hi G H( S+2 \X) and obtain an 
estimate for 

\\d z >G( z (2), z} (u ) - d z 'G( z (i), z )(u a )\\ H{s) 
of the same form as in l|3.17|) . □ 

The two previous lemmas yield 

Proposition 3.5. Let s G R, a subdivision of [0, Z] as in Definition lff.il 
and let u a G H^+^iX). Then the map W«p, z (uo) is ^f°([0, Z], and 
piecewise ^([O, Z], H^>(X)) is chosen such that Atp is small enough. The 
map z 1— ► Wfp )2 (iio) is in fact globally Lipschitz with C > such that 

||W«p, z '(u ) - yV<p, 2 (uo)||,H(=) < C\z' - z\\\u \\ H{ s+i). 

We recall that U(z', z) is the solution operator of the Cauchy problem Q1.5|) - 
(|1.6() . We can then apply the energy estimate (|1.7|) to U(z,0)(uo) — W<p,z(uo) 
(adapt the proof of Lemma 23.1.1 in [Hj to the case of a Lipschitz piecewise C 1 
function) and obtain 

(3.18) sup exp[-Az] ||f/(z, 0)(w ) - V%.*(u )||jf<.> 

z€[0,Z] 

<2 [ exp[-Az] \\(d z + a z (x,D x ))yV VtZ (u )\\ H( s)dz. 
Jo 

Let M G H( S+1 1(X) and let qj = {z<°\ . . . , z^}. We take z e]z (fc) ,« (fc+1) [. 
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Then 

(d z +a z (x,D x )) V%, z (u ) 

= (d z + a z {x,D x )) {^( z , z w)\\Q( z ^). z ^-i)){u ^ 
= (d z + a z {x,D x )) (Q( z , z m)(uk)) 

k 

with u k := Y[G( z (i), z a-i))(u ) which is in by Theorem E22 We 

i=l 

first turn our attention towards the term (d z + a z (x, D x )) \Qi z z (k)\(u)j for any 
u G H^ S+1 \X) as the norm of Uk in H^ S+1,> (X) remains under control even if 
|5)J| = N becomes very large by Proposition 13. 21 

(3.19) 3K>0, \\u k \\ H(s+1) <K\\u \\ H(s+1) , ke{0,...,N}, 

n= |«p| gN, u e ff (s+1) pc), 

if Asp is small enough. 

We shall need the following lemma which is a variant to Lemma fe/ill 

Lemma 3.6. Let q A (z,x,y,£) be an amplitude in Sz(SL p x W) depending on 
the parameters A > and z G [0, Z] that satisfies Property \Ql\I o,nd such that 
Qa(z,.)\a=o = 0. Let r(x,y,Q G S s (R 2p xRf) for some s 6l. T/ien 

(7{gA rj^^O-^A^^^O r(a:,a:,0 = A m+2i A A l (z, x, 0, < m < p - 5, 

where the function X l A (z,x, £) is bounded with respect to A and z w»t/i values in 

Proof. We proceed as in the proof of Lemma 12.211 (we take p — 1 for the sake 
of concision) . We obtain 

cr{q A r} (z,x,£) - q A {z,x,x,£)r(x,x,X) 

= ° I ® 3 ® 4 ^ A ( Z ' x ' ^ — s ^ x ^ ^ s 
where here 

939 4 A A (2:,a;,y,C) = (d y d^q A )(z,x,y,£,) r(x,y,£) + d y q A {z,x,y,t)d e r(x,y,t) 
+ d i q A (z,x,y,£,)d y r(x,y,^) + q A (z,x,y,£,)d y d s r(x,y,£,). 

The first two terms are treated like in the proof of Lemma [2.211 For the Third 
term, with Property \Ql\ we write 

%A d y r = A m ' +S g™' (00)1 8 y r, < mf < 1 - S 

where q^ { ° m d y r G S , ™' +s -' , (]R 2 p x W). We actually take <5 < to' < 1 - 6 and 
write to = to' — 5. We obtain 

£l „ .0 ™ A rn+25 ~m 

d^q A O y r = A g A , 
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where q% is bounded w.r.t. A with values in s™+ s -p+ s (R 2 p x W) and < m < 
1 — 28 = p — 8. For the fourth term we write 

Am' m' (0,0)0 n ^ I ^ i 

g A = A q A < m < 1, 

where g ™' (0 ' 0)0 e S™'(R 2 f> x MP) by Property (glj since <? a |a=o = 0. We 
actually take 28 < ml < 1 and write m — m! — 28. Then 

q A d y d s r = A m+2S q A \ 

where q A is bounded w.r.t. A with values in 5™ +i ' _ ' p_ '^(R 2p x W) as m+ s — 
1 + 28 = m + s — (p — S) and < to < 1 — 28 = p — 8. We conclude like in the 
proof of Lemma \2. 2 II □ 

For the next proposition we shall need the following assumption as announced 
in Section ^ 

Assumption 3.7. The symbol a(z,.) is assumed to be in ([0, Z], S 1 (M. n x 
R n )) } i.e. Lipschitz continuous w.r.t. z with values in S 1 (M n xK"), in the sense 
that, 

a(z', x, £) — a(z, x, £) = (z! — z)a(z' , z, x, £), < z < z' < Z 
with a(z',z,x,£) bounded w.r.t. z' and z with values in 5' 1 (R n x R ra ). 

Proposition 3.8. Let s G R. There exists Aq > and C > such that for 
z' -z = A, Ae [0,A 6 ], 

+a z '(a;,-Dx))0( z ', 2 )||(iTW,/f(.-i)) < CA^ 

Like in the proof of Theorem 12.221 we assume that c\ satisfies property Pl for 
some L > 2. We know that it is always true for L = 2 by Lemma [2.131 but 
special choices for c\ can be made. As before we use p = 1 — 1/L and 8 = 1/L 
with p > 8 for L > 2 and p = 8 = \ for L = 2. 

Proof. Let At z i tZ -\ be d z >G( z < yZ ) an d ^(z',z) be a z i(x,D x ) oQ^ z , z ^ with respective 
kernels A^ z i^(x' ,x) and B^ z i z ^(x',x). We have 

Az'.z) ( x 'i x ) = ~ J exp[i(ar' - 2c|£)] exp[- Aa(z, a/, £)] a(z, or', £) oT£. 
Let us define 

£>(z',z) := (-4(z', z ) + ° ^ _2s ° (-4(2', «) + ,*))*• 

We prove in the following lemma that for r, s S R, ||2'(2',2)||(ffW,ij(r+2s-2))) < 
CA uniformly w.r.t. z G [0, Z] for A small enough. The conclusion then follows: 
if C (zV) := £ s_1 o Z>(a/ )Z ) o i; 8 - 1 then ||C (zV) < CA (take r = -s+ 1); 
then H^- 1 o (^ (z ^ } + 6 (2V) ) o ^- s || (i2!i2) < CAi □ 

Lemma 3.9. Let r,s G R. Then W'D^'^Wih^ ,fl"( r + 2s - 3 ))) — uniformly 
w.r.t. z G [0, Z] /or A small enough. 
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Proof. The operator T>/ z / z \ is made up of four terms: 

:= -A {z ,, z) o £T 2s o A* (z , tZ) , Vu*>,z) ■= Az\z) o E- 2s o B\ z ,^ 
V 3 , ( z>,z) ■= B {z ,. z) o E- 2s o P 4i(2V) := o E- 2s o B^. 

The kernel of fi,^'^) is given by 

= J exp[i(x' - x\§ +iA{b 1 (z,x',C) ~ h{z,x,£))] di, z (x' ,x,£) <2£ 

where 

di,z0',2:,£) = w (z / iZ) (a;',ar,^) a(z,a;',£) 

and 

a;, := g {z ^ z) (x, f) g^ z >, z )(x, £) exp[-A(ci(z, a/, £) + a(z, x, £,))} 

with g( z ',z) given in (|2.9[) . Following the proof of Theorem 12.221 we write 
bi(z,x',£) — bi(z,x,£) = (x 1 — x\h(z, x', x, £)) where h is homogeneous of de- 
gree one in £, |£| > 1. The function h and continuous w.r.t. z with values in 
S 1 (X x R n ). We thus obtain that the change of variables £ — > £ + A/i(2, a;', a;, £) 
is a global diffeomorphism for A small enough (uniformly in z <E [0,Z]). The 
Jacobian J7a(z, x', x, £) is homogeneous of degree zero in £, w.r.t. A and 
bounded w.r.t. z with values in 5°(R 2n x M. n ). We then have 

Di, {z ,, z) {x',x) = J exp[i(a/-a:|£)]di, z (^z,£(A,£)) J(A, z,x' ,x,£) <f£. 

The function £(A, 0, x', x, £), written £(A,£) for concision, is bounded w.r.t. z 
and < ^°° w.r.t. A in S' 1 (R 2 ™ x M") and homogeneous of degree 1 in £ as shown in 
Lemma 12.231 It follows that di yZ (x\ x, £(A, £)) J(A, z, x' , x, £) is then bounded 
w.r.t. z and A with values in S 2 ~ 2s (W 2n x R") by Lemma 12.101 and the proof of 
Theorem l2~2~2l Note that if A = then f(A,£) = £. The operator X>i,( z ', z ) is 
thus in vl/ 2_2s with symbol 

rfi,( 2 ',z)(a;',C) = o-{di, z (x',x,£(A,£)) J{A,z,x',x,Cl) 
Similarly we prove that A( z ' yZ ) o E~ 2s o Q* z , s is the ij)T)0 with amplitude 

-u {z , tZ) (x',x,i(A,$) a(z,x',i(A,0) J(A,z,x' ,x,£). 
The operator £>2,(z',z) is thus in ^> 2 ~ 2s {X) with symbol 

d>2,(z' ,z)(x' , 

= -a {up,,) (a/, a, f(A, 0) a(z, a/, f (A, £))J(A, z, x', x,{)} # a*(z', x', £) 
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Similarly we find that the operators T>s,(z',z) an d T^a,{z',z) are i n 2s (^0 with 
respective symbols 

d 3 ,(z', z )(x' ',0 = -a(z',x',^) 
and 

d4,(*',,)(a / .0 = o (« / . aj/ .0 J(A,« 1 a/,a: J e)} 

#a*(z',x',0, 

For g(a;', i, £) an amplitude we define 

Z{q}(x',0 :-a{(e>- 2s a{z,x>,Q q(x',x,0 a(z,x,0) 

-a{(0- 2s a(z,x',0 q(x',x,0} #a*(z,x',Q 
+ a(z,x',0 # a{(0- 2s q(x'x,H)U a*(z,x',0 

-a(z,x',0#°{{t,r 2s Q{x'x,£) a(z,x,£)}. 
The operator T>^ z / ^ is thus in ^> 2 p ~ 2s (X) with symbol 

d(z',z) — dl,{z',z) + dn,(z',z) + ^3,(z',z) + di,(z',z)- 

Such a symbol is bounded w.r.t. A, for A small enough, as the composition 
formula for symbols is a bounded map. Note that 

9(z',z){x',0 3(^(x,£)(0" 2s i7(A,z,x',:e,0 = (0~ 2s + Ak A (z,x',x,£) 

with kA bounded w.r.t. z and < ^ , °° w.r.t. A with values in S~ 2s (X' x X x R n ) 
as oj/ z t fZ \J is itself < ^°° w.r.t. A and 9( Z ',z) <7(z',z)»7|a=o = 1- By Assump- 
tion we also write a(z',a;,£) = a(z,x,£) + Ao(z',z,2;,() with a(z',z,.) 
bounded w.r.t. z and A with values in S 1 (X x R). We thus obtain 

V(z',z) =V a (z , iZ) + AVl, >z) 

with symbols 

d? z , tZ) :=X{{0- 2s Pa(z,x',x,0}, 



and dh, . which is bounded w.r.t. z and A with values in S 2 2s (X' x X x R"). 
The Calderon-Vaillancourt theorem (see ^| Chapter 7, Sections 1,2] or 
Section XIII-2]) in the case L = 2 or Theorem 18.1.11 in [Sj in the case L > 2 
yields \\T>} Z , z - ) \\(H( r \H( r + 2 °- 2 ))) < ^l- Note that for a symbol q(x',£) we have 
£{<?(x',£)} = as 

a{<z(x',0 n(z,x,£)} = q(x',£) # a*(z, x', = a{q(x', £)} # a*(2,x',£), 

for any symbol q. Thus d a (z',z) = Y^{(^)~ 2s (pa(z, x', x, £) — 1)}. Lemma [3.61 
allows us to write (take m = p — 5) 

ff{(0~ 2s (PA(2,a;',2;,£) - l)a(z, x', £)a(z, a?, 0} 

= (0~ 2s Cpa(z, x', x', - l)o(«, x', £)3(z, x', £) + AA A ,i(«, x' , f) 
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where Aa i bounded w.r. t. z and A with values in S 2 ~ 2s {X' x R"). We also 
write 

<t{(0- 2s (pa(z,x\x,0 - l)a(z,x',0} # a*(z,x',0 

= ((0- 2s (pa(z,x',x',0 - l)a(z,x',0) # a*(z,x',0 
+ AX A:2 (z, x', £) # a*(z, x', £) 

= <?{ (0 ~ 2s (PA (z, x',x',0- x', £)a(z, x, £)} 

+ AX A<2 (z,x',0 #a*(z,x',S) 
= (0~ 2s (pa(z,x',x',£,) - l)a(z,x',£,)a(z,x',£) 

+ A(X A>3 (z, x', + Aa )2 (z, a/, C) # a* (z, x', £)) 

where Aa,2 and Aa,3 are bounded w.r.t. z and A with values in S^~ 2s (X' x R") 
and S 2 ~ 2s (X' x R n ) respectively. Similarly we have 

v{(£,y 2s (PA(z,x',x,£,) - 1)} # a*(z,x',£) 

= ((e>- 2s (PA(^.x'^',0 - 1)) # a*(«,a/,0 + AAa,4(^x',0 # 
= a{<0" 2s (PA(^, a;', x', £) - l)a(z, x, £)} + AA A , 4 (z, x', f ) # a*(z, x', £) 
= (i)- 2s { P A(z,x\x' 1 0-l)a{z 1 x',0 

+ A(A A , 5 (z, x', + A Al 4(z, a/, # o*(2, x', £)) 

where Aa.4 and Aa,5 are bounded w.r.t. z and A with values in S~ 2s (X' x R") 
and 5'p _2s (X' x R") respectively and 

<7{ (0 - 2s (PA 0-1)^,^,0} 

= (0- 2s (pa(z,x',x',0 - l)a(z,x',C) +AX A , 6 (z,x',0 

where Aa.6 bounded w.r.t. z and A with values in S p ^ 2s (X l x R"). We thus 
obtain 

d a {z , z) = A(A A ,i + A A , 3 + A A ,2 # a* + a # A Aj5 + a # A A . 4 # a*) = 

with , z s bounded w.r.t. z and A with values in S 2 ~ 2s (X' x 1"). This con- 
cludes the proof. □ 

We have thus obtained a convergence result in the Sobolev space fTW(M n ) for 
Wqj iZ (uo) if the initial data uq is in ij( s+1 )(R n ). The result is actually the 
convergence of the Ansatz Wsp, z to the solution operator U(z, 0) in the norm of 
L(#( s+1 )(R n ),.ff( s )(R n )): 

Theorem 3.10. Assume that a(z, .) is in if ([0, Z], S 1 (M n xR n )), i.e. Lipschitz 
continuous w.r.t. z with values in S 1 (W n x R"), in t/ie sense that, 

a(z , x, £) — a(z, £) = (z' — z)a(z', z, x, £), < z < z < Z 

with a{z' , z , x , £,) bounded w.r.t. z' and z with values in S' 1 (R n xl"). Let s E R. 
TTien the approximation Ansatz W«p iZ converges to the solution operator U(z, 0) 
of the Cauchy problem £^-£1$ in L(ff( s+1 ) as A<p 5 oes to 

with a convergence rate of order ^ : 

\\W^, z -U{z,Q)\\{H^),HM) <CA|. 
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Proof. Using (|3.18[) and l|3.19[) we obtain 

sup exp[-Az] \\U(z,0)(u ) -W<p lZ (u )\\ H (,s) 
ze[o,z] 

< 2 f exp[-Az] AlCK\\u \\ H(s+1) dz < Ca|||uo||jj<.+i). 
Jo 

The result follows. □ 

If we change the assumption made on the symbol a(z, .) to some Holder type 
continuity, then the corresponding change in the proof of Lemma y . 91 yields the 
following weaker result 

Theorem 3.11. Assume thata(z,.) is in ^ a ([0, Z], S 1 (W l xl")), i.e. Holder 
continuous w.r.t. z with values in S 1 (M n x R n ), in the sense that, for some 
< a < 1 

a(z , x, £) — a(z, x, £) = (z — z) a a(z , z, x, £), < z < z < Z 

with a(z' , z, x, £) bounded w.r.t. z' and z with values in S 1 (W l x M"). Let s 6 M. 
T/ien i/«e approximation Ansatz W<p jZ converges to the solution operator U{z, 0) 
o/tte Cauchy problem EJft-fTQil in L(i7< s+1 ) (M™), ffW(R n )) as A<p 5 oes to 
wif/i a convergence rate of order a/2: 

\\W v , z -U(z,0)\\(h(^Kh^)<CA« /2 . 

A result similar to that of the previous theorems can be obtained with weaker 
assumptions, namely without assumptions on the symbol a(z, .) like those made 
in Theorems 13.101 and 13.111 by introducing another, yet natural, Ansatz to 
approximate the exact solution to the Cauchy problem (|1.5fl - l|1.6fl . For a symbol 
q(z,y,r)) € tf°([0, Z], S m (RP x R r )) we define q {z >, z) (y, n) £ <g°([0, Z} 2 , S m (RP x 
R r )) 



%z',z)(y,v) ■= — / q{s,y,ri) ds. 

z z J z 

Then we define 

(3.20) ^, g) (x',x,0 :=(x'-x\0+iAa 1{zf}Z) {x',0 

= (x'-x\£) + Ab 1{z , iZ) (x'^) +^c 1{z ^ z) {x',0- 

and 

(3-21) 9(z',z){x,Cj := exp[-Aa ( z ' iZ )(x,^)}. 

and finally, following [T7|. we denote G( z ',z) the FIO with distribution kernel 



G^'.z) (x', x) = J exp[i(x / -x|^)]exp[-Aa( z / )Z )(a; / ,^)]3^ 

= / ejq>[i$( z , jZ )(x',x,£)}g( z , tZ) (x',£)fi£. 



with the associated approximation Ansatz 
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Definition 3.12. Let = {z^\ z (1 \ . . . , z^} be a subdivision of [0, Z] with 
= < < • • • < = Z such that z^ +1 ~> - = A v . The operator 

Wfp, z is defined as 



W<p, 2 := < 



Q{zfi) if 

k 

W))!!^),.!*- 1 )) # zik) ^ z ^ z(k+1) - 



Most results of Sections and EJapply to this new Ansatz. We give some details 
about how to adapt some of the proofs. We have 

Lemma 3.13. Let q{z,y,rj) £ ^ ([0, Z],S 1 (W x K r )) that satisfies Property 
(Ftp. Then q< z ' jZ \{y,rj) also satisfies Property \Pl\) . 



Property \Pl\ in Definition 12. Ill is now to be understood w.r.t. to two parame- 
ters z' and z. 

Proof. Uniform bounds w.r.t. z and z' will be immediate. The case \a\ + 1/3| > L 
is clear by Remark 12.121 Let then \a\ + \(3\ < L and observe that 



\dX^'Mv,v)\ = l-^r- / dytfqisMfi) ds\ 

z z J z 

<(7(l + |r7|)^l+(H+l/3|)/£_L_ r (l + g(z, 2 ;,ry)) 1 -(l Q l + l' 3 l)/ i d S 



/ < v l-<|«| + |/9|)/i 

<C(l + |r;|)-l' 3 l + (l«l + l' 3 l)/ L 1 + ^— / q(z,y,v)ds) 



= (7(1 + I^D-l^l+d-l+I^D/^l + ^ 7?)) l-(|a|+l^l)/i 

by Jensen inequality as t —(l+t) 1_ (l Q l + l /3 l^ L is convex when + < L. □ 
As a consequence of Lemma 12.171 we have 

Lemma 3.14. Let q(z,y,rj) £ ^([O, Z], S 1 (W x W)) that satisfies Property 
Then p A := exp[— Aq^ z ,^{y, Tj)] satisfies Property jQ L \) - 



The result of Theorem 12.221 thus applies to the modified thin-slab propagator 
G( z ',z) (Lemma |2.23l has to be slightly modified). The proof of Lemma l3.3l applies 
with the aid of Proposition 12 .251 as 

(G(zV),z) ~ G(zW,z)){uo){x') = 
Z<2) ff 

// exp[i(x —x\Q — (z — z)a/ z i >z \(x ,£)] a(z ,x ,£) Uo(x) dx cf£ dz . 

(!) J J 

To adapt the proof of Lemma 13.41 we need 

Lemma 3.15. Let s € M and z" , z £ [0, Z\. The map z' i— ► d z iQt z i <z \, f or 
z' £ \z" ,z\, is continuous with values in L{H^ S+2 \X),H^ S \X)), forz" — z = A 
small enough. 
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Proof. We choose A = z" — z sufficiently small such that the results of Section^] 
apply. Let z^ x \ z^ 2 ' G [z,z"]. Then we have 

d z >G (z (2) tZ) (x',x) - d z /G( z (i) jZ) (x',x) 

= - / ~ -101 (^ (2) . o «p[-/f «(•. *] 

-a(z {1 \x\ £) exp[-J z z a(s, x' , £) ds] \ d£ 

= A (z (2) z (i) z )(x',x) + B( zi 2) z (y ) z j(x', x), 

where 



A (zV\z(»,z)i x ', x ) ■= - J exp[i(x' - x\£)] a(z (2) ,x' ',£) 

/ z (2) z (l) \ 

(exp[-/ z a(s,x',£) ds] - exp[-/ a a(s, x', £) ds}) d£, 

and 

{a{zW,x>,Z)-a(zW,x>,S)) exp[-/f ^(s, x', ds] d£. 

We write 



A( z (2) iZ m >z) (x',x) = J ^ J exp[i{x' - x\£)] 

a(z (2) ,x',£) a(s,x',£) exp[-(s - z)a( s<z )(x' ,£)] ds d£. 

and for the associated operator, ^4( z (2> z (i> z \ we obtain by Proposition 12.251 
that ||yl( z (2) , z (i). z ) < C|z( 2 ) — For the second term we can 

apply Proposition 12.251 which gives the estimate, for the associated operator, 
||2?( z (2) jZ (i) jZ )(a/, x)\\( H (s+2) H (s)j < C p(a(z( 2 \ .) — a(z^\ .)), with p a seminorm 
in ^(A x R n ). The continuity of z ^ a(z, .) in S 1 (X x M n ) (Assumption EJ 
yields the result. □ 



With the previous lemma we can easily adapt the proof of Lemma l3 . 41 and obtain 
the same result for G( z ', z )- 

Lemma 3.16. Let s G E, z", z G [0, Z], wii/i z < z", and let u G # (s+1) (X). 
TTiera ifte map z' i-> ^, z )(tt ) is in c ^([z ) z"] 1 H^+ 1 \X))C\^ l {[z,z"] 1 H^\X)) 
for z" — z = A small enough. 

This allows to use the energy estimate (|1.7ll . 

We now note that in the proof of Lemma 13.91 with the new thin-slab propa- 
gator, G( z '. z ), the amplitudes of the operators T>\, . . . ,2?4 only involve the term 
a(z', x, £) instead of both a(z', x, £) and a(z, x,£) (as d z /((z' — z)ai z r iZ \(x' , £)) = 
a(z', :£,£)). Thus the proof of Lemma 13.91 does not require any assumption like 
Assumption 13.71 made in Theorem l3.1Ul or assumptions of Holder type regularity 
on the symbol a(z, .) made in Theorem 13. Ill Consequently we obtain 
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Theorem 3.17. Let s£l. Then the approximation Ansatz W«p, z converges in 
L(H ( - a+1 '>{W l ),HW{W 1 )) to the solution operator U(z,0) of the Cauchy prob- 
lem fT5)) - f71>)) as Aqj goes to with a convergence rate of order i : 

\\W v . z -U(z,0)\\ {His+1)Mis)) <CAl. 

A A diagonalization/decoupling of the acoustic 
wave equation 

We give here an overview of |22| . which gives a motivation for approximat- 
ing solutions of the Cauchy problem l|1.5|) - (|1.6|) . f° r instance in the context of 
geophysics. 

We first consider the scalar wave equation 

(A.22) (^-p-'c-'d? + J2 djp-'d^j u = F, 

as encountered in acoustics, where p is the fluid density, and c is the wavespeed. 
Both these functions are assumed to be independent of time t and to be in 
< ^°°(R"). We further assume that < p < p{y) < Pi and < c < 
c (y) — c ii V ^ We denote z = y n and x = (j/i, . . . , y n -i) and write 

p(x,z,D t ,D x ,D z ) = p^c-^Dl - YT j ZlD j p~ 1 D j - D z p~ x D z where D = \d. 
Its principal symbol is P2{t, %, z, r, £, C) = j o _1 (c _2 r 2 — |£| 2 — Q 2 ). 

Note that r ^ in Char(p). We put lfA~22)) in a matrix form 

(A.23) D z w(t, x, z) = G(x, z, D u D x )w(t, x, z) + f(t, x, z) mod < ^°°, 

«*°-{°a A o')' -(AO- '-{})■ 

where A is a first-order elliptic ipDO, say for instance |Z) ijSC |, and 

n-1 

A = p^c^DlK- 1 - DjP^DjA' 1 , 

with A -1 denoting a parametrix for A. 
Following |22] j we introduce 

I' @ = {(x,z,t,0 I t^0,\c(x,z)t-^\ <sin9}, 

Ie = {(t,x,z,r,£,,() I (x,z,t,0 e I'e, ICI < c^Vlh 

where O € (0, -|). The inequality |£| < c(x, z) _1 |r| on Char(p) explains the 
condition |£| < c^lrl above. We choose an angle O 6 (0, ^) and work in the 
microlocal region Iq assuming that WF(w) C Iq. FigureQilhistrates the set Iq 
at a given (x, z) and a given frequency r. An angle 9 £ [— O, 0] corresponds to 
a propagation angle. Restricting the analysis to Iq corresponds to staying away 
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Figure 1: The shaded area corresponds to Iq at a given (t,x,z) and a given 
frequency r. is the propagation angle. The set Char(p) is represented dotted. 



from horizontal propagation. Note that in Iq we have c(x,z) 2 r 2 — |£| 2 > 0, 
which is the main purpose of the restriction to such a microlocal region. 

In Iq, G is a first-order ^DO by Theorem 18.1.35 in 8 . In Iq we can follow the 
method of [23 Chapter IX] (see also |24|) to decouple the up-going and down- 
going wavefields. We briefly recall the method here. Define rf il {x, z, r, £) = 
±(c(x, z)~ 2 t 2 — |£| 2 )2, which are the two roots of det(7y/2 — Gi) = with G\ 
the (matrix-) principal symbol of G. The matrix G\ (x, z, r, £) is diagonalizable 
and we choose a matrix V(x,z,t,£) E S°(I' e ), invertible, such that VGiV^ 1 
is diagonal; V van be chosen homogeneous of degree 0. If we write = 
V(x, z, D t , D x )w we obtain 

,x,z) , x, z) mod ff 00 , 

G (0) = (D Z V)V- X + VGV' 1 mod in Iq, /<°) = V f. 

We write G<W = G ( ° ] + g[, 0) with G^ 0) G tf 1 in J e and diagonal and g[, 0) G *° in 
Iq. By V^ _1 we denote a parametrix for V with principal symbol V(x, z, r, £) _1 
(an abuse of notations, which will occur below again). 

We then write u;W = (1 + K^(x, z, D t , D x ))w^ , with G f" 1 in Iq of 



the form 
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We then obtain 

D z w^ = G^w^ + [KV>,GW]wW + G^wW + + R w w^ mod <^°°, 

g *-i in/ 0) /« = (l + K( 1 ))/(°), 

making use of 

(1 + K^)Gf ] (1 + iff 1 ')- 1 = G[ a) + [ifW, G< 0) ](1 + K W)- 1 

and the fact that L(l + Iff 1 ))- 1 - L € S^" 1 if X € * m Lemma 2.1 in [21 
shows that if W can be chosen so as to have [K^\ G^] + Gq ' diagonal up to 
an operator in 'f -1 in I@. The procedure goes on by choosing S in 
7e in order to diagonalise the term of order -1, etc. We thus obtain Q S ^° in 
7e such that w = Q^ 1 w satisfies 

D z w^Gw + f mod-r 30 , f^Q- 1 /, 
with G = G(x, z, Z?t, 1)^) € , J /0 in ie, diagonal up to a regularizing operator 

G - ^ 6_ 

In [221, Stolk shows that b± can be chosen selfadjoint. This is achieved by first 
choosing selfadjoint operators with principal symbols equal to Ty ± (x, z, t, £) and 
then replace (1 + K^) by exp[if M] in the iteration process described above. 
Various choices of Q are presented in |22|. 

We define the set Je+ of points (to, xq, zq, t , £o, Co) such that the bicharac- 
teristics associated with b+, parametrized by z, (t(z),x(z),T(z),£(z)), pass- 
ing through (to, Xo, To, £o) at z — zq, is such that for all z € [0, Z], the point 
(x(z), z, £(z), t(z)) remains in I'q. In other words, with the interpretation given 
by Figurenjthe propagation angle, 9(z) along the bicharacteristics should never 
exceed 6. 

We now choose < ®\ < 02 < f . We choose a real non-negative symbol 
c(z,x,t,£) € S 1 (R x R™- 1 x K x R"" 1 ) such that c = in J 01 and elliptic in 
the complement of I® 2 . After extendig smoothly b+ outside Is , such that b + is 
real homogeneous of degree 1, we now consider the Cauchy problem 

(d z - ib + (z, x, D t ,D x ) + c(z, x, D t ,D x ))v = 0, 

v(0,.) = v+(0,.), 

where 

With Assumption (33) and (34) in |22| we obtain that 

v = v + mod "if 00 in Jei+j 

f = mod ( ^' 00 in the complement of Je 2 +- 

See [22 and [2U for details. A similar results holds for the other 'one-way' wave 
operator d z — ib- + c. 
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